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ABSTRACT 

Context. Starless cores represent a very early stage of the star formation process, before collapse results in the formation of a central 
protostar or a multiple system of protostars. 

Aims. We use spectral line observations of a sample of cold dust cores, previously detected with the BLAST telescope in the Vela-D 
molecular cloud, to perform a more accurate physical and kinematical analysis of the sources. 

Methods. We present a 3-mm and 1.3-cm survey conducted with the Mopra 22-m and Parkes 64-m radio telescopes of a sample of 
40 cold dust cores, including both starless and proto-stellar sources. 20 objects were also mapped using molecular tracers of dense 
gas. To trace the dense gas we used the molecular species NH3, N^PP, HNC, HCO + , H 13 CO + , HCN and H 13 CN, where some of them 
trace the more quiescent gas, while others are sensitive to more dynamical processes. 

Results. The selected cores have a wide variety of morphological types and also show physical and chemical variations, which may 
be associated to different evolutionary phases. We find evidence of systematic motions in both starless and proto-stellar cores and 
we detect line wings in many of the proto-stellar cores. Our observations probe linear distances in the sources J 0.1 pc, and are thus 
sensitive mainly to molecular gas in the envelope of the cores. In this region we do find that, for example, the radial profile of the 
N 2 H + (1 - 0) emission falls off more quickly than that of C-bearing molecules such as HNC(1 - 0), HCO + (l - 0) and HCN(1 - 0). We 
also analyze the correlation between several physical and chemical parameters and the dynamics of the cores. 

Conclusions. Depending on the assumptions made to estimate the virial mass, we find that many starless cores have masses below 
the self-gravitating threshold, whereas most of the proto-stellar cores have masses which are near or above the self-gravitating critical 
value. An analysis of the median properties of the starless and proto-stellar cores suggests that the transition from the pre- to the 
proto-stellar phase is relatively fast, leaving the core envelopes with almost unchanged physical parameters. 

Key words, submillimeter: ISM — stars: formation — ISM: clouds — ISM: molecules — radio lines: ISM 



. £h ■ "'" ' n t r OduCtion rived luminosities and masses. Recent surveys with the MIPS 

k> . instrument of the Spitzer Space Tele scope are able to c onstrain 

j_j - Starless cores represent a very early stage of the star formation the temperatures of warmer objects (ICarev et al.ll2005l) . but the 

d [ process, before collapse results in the formation of a central pro- youngest and coldest objects are potentially not detected, even 

tostar or a multiple system of protostars. The physical proper- in the long-wavelength Spitzer bands. 

ties of these cores can reveal important clues about their nature. Qn ^ Qther hand ^ more recent g wkh both 

Mass spatial distributions, and lifetime are important diagnos- ^ BLAS T ("Balloon- borne Large-Aperture Submillimeter 

tics of the main physical processes leading to the formation of Tdescope » p ascale et al 2008) and Herschel t elescopes (see, 



the cores from the parent molecular cloud. Netterneld et al. 2009. 1oimi et al.l 120091 iMolinari et alJ 



In the past, (sub)millimeter continuum surveys performed 1201 Ol) have demonstrated the ability to detect and characterize 

with ground-based instruments have probed the Rayleigh- Jeans cold dust emission from both starless and proto-stellar sources, 

tail of the spectral energy distribution (SED) of these cold ob- constraining the low temperatures of these objects (T Z> 25 K) 

jects, far from its peak and, at short submillimeter wavelengths, using their multi-band photometry near the peak of the cold 

have been affected by low sensitivity due to atmospheric condi- core SED. The bolometric observations alone, however, cannot 

tions. Therefore, these surveys have been limited by their sensi- fully contrain the dynamical and evolutionary status of the cores, 

tivity or by their relativ e inadequacy to measure the temperature Spectral lines follow-ups are necessary in order to investigate the 

(e.g jMotte et alJ|j998h . producing large uncertainties in the de- physical, dynamical and chemical status of each detected core. 
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Table 2. Molecular species observed with the Parkes and Mopra 
telescopes 



1 - 



Spectral line 


E u 




Rest Frequency 




[K] 


[xlO 3 citT 3 ] 


[MHz] 


NH 3 (1, 1) 


23.4 




23694.496 


NH 3 (2,2) 


64.9 




23722.633 


N 2 H+(1 21 - 32 ) 


4.47 


2.5 


93172.053 


HNC(l 2 -0i) 


4.35 


9.5 


90663.574 


HCO + (l -0) 


4.28 


2.5 


89188.526 


H 13 CO + (l - 0) 


4.16 


2.2 


86754.330 


HCN(l 2 -0i) 


4.25 


28 


88631.847 


H 13 CN(1 2 -0i) 


4.14 


20 


86340.167 



Notes. (fl) See Section l4~4~2l for a definition of n eS . 



BLAST has identified more than a thousand new star- 
less and proto-stellar cores during its second long duration 
balloon science flight in 2006. BLAST detected these cold 
cores in a ~ 50 deg 2 map o f the Vela M olecular Ridge 
(VMR) (INetterfield et al I 120091) . The VMR dMurphv & Mavl 
1991, iLiseau et al.lll992l) is a giant molecular cloud complex 
within the galactic plane, in the area 260° ^ / 5= 272° and 
-2° 5: b $ 3°, hence located outside the solar circle. Its 
main p rop erties have been rec ently revisited by INetterfiel d et alj 
(120091) and lOlmi et alJ (l2009h 

In particular, lOlmi et alJ (|2009) used both BLAST and 
archival data to determine the SEDs and the physical parameters 
of each source detected by BLAST in the smaller region of Vela- 
D (se e Figure Q]), where observati ons from IR dGiannini et all 
2007) to millim eter wavelengths dMassi et al. 2007) were al- 



ready available. 



data of Eliaetal 



Olmi et all d2010h then used the spectral line 
(120071) analysis of the dynam- 

ical state of the cores toward Vela-D. They found that ~ 30-50% 
of the BLAST cores are gravitationally bound, and that a signifi- 
cant number of cores would need an additional source of external 
confinement. 

In this work we will use spectral line observations, both 
single-points and maps, of a sample of BLAST sources from 
Vela-D to perform a more accurate physical and kinematical 
analysis of the starless and proto-stellar cores. This paper is thus 
organized as follows: in Section[2]we describe the selected tar- 
gets and molecular lines. In Section[3]we present the results of 
the analysis of the observed molecular line spectra and maps, 
describing the derived physical and kinematical parameters. We 
then further discuss these results in Section |4] and conclude in 
Section|5] 



2. Observations and archival data 

2.1. Targets and lines selection 

The observations were performed in June 2009 with the Mopra 
and Parkes telescopes of the Australia Telescope National 
Facility (ATNF). Althou gh original l y our aim was to observe all 
sources in the catalog of Olm i et aD (l2009l) . because of time con- 
straints we were able to observe a total of 40 sources at Mopra 
and 22 sources at Parkes, in either mapping or single-pointing 
mode (or both, see Table Q] and Figure [TJ. We have selected an 
almost equal numbe r of starless (19 ) and proto-stellar (21) cores 
where, following Ol mi et al. (2009), we define the BLAST cores 
as proto-stellar when they turned out to be positionally asso- 
ciated with one or more MIPS sources at 24-fim, and starless 
otherwise. Each sub-sample is composed by objects that present 
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Fig. 1. The gray-scale image shows the BLAST 250/im map 
of Vela-D, with galactic coordinates in degrees. Superimposed 
are the locations of both starless (ope n circles) and proto-stellar 
(crosses) cores (see Ol mi et al. 2009). The size of the symbols 
indicates the mass range [M ] of each core and color-coding in- 
dicates the core temperature [K] (see legend). 



increasing temperature and different morphology. This selection 
criterion is quite arbitrary but aimed at identifying different evo- 
lutionary stages within each sub-sample. 

In order to probe the dense and cold gas detected by BLAST 
toward the cores in Vela-D we needed appropriate molecular 
tracers. We thus selected several rotational transitions with low 
quantum numbers to be observed at Mopra in the 3-mm band, 
and we also observed the NH3(1,1) and (2,2) inversion lines with 
the Parkes telescope. In the following sections, our analysis will 
be based on the ammonia lines and on the main spectral lines ob- 
served at Mopra, i.e., N 2 H + (1 -0), HCN(1 -0), HNC(1 -0) and 
HCO + (l - 0) (and two of their isotopologues, see Table|2]i. This 
choice of molecular tracers ensured that we would be able to 
trace the dense as well as the more diffuse gas, and we would 
also be able to detect molecular outflows. Furthermore, these 
molecules should also be able to give out clues about the dif- 
ferent chemical states of the sources. 



2.2. Mopra 

During the spectral line mapping or single-pointing mode ob- 
servations with the Mopra telescope the system temperatures 
were typically comprised in the range ~ 170 - 280 K and we 
reached a RMS sensitivity in T* units of about ~ 30-40 mK 
after rebinning to a 0.22 kms~' velocity resolution. The beam 
full width at half maximum (FWHM) was about 38 arcsec and 
the pointing was checked every hour by using an SiO maser as 
a reference source. Typically, pointing errors were found to be 
~ 5 - 10 arcsec. 
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Table 1. BLAST sources in Vela-D observed with the Mopra and Parkes telescopes 





BLAST 






Mopra 


Parkes 


Source W 


Source name 


/ 


b 


Starless 


Observation 


Observation 






[deg] 


[deg] 


or Proto-Stellar 


type* 


type 6 


3 


BLAST J08453 1-435006 


263.6001 


-0.5302 


P 


SP 


SP 


9 


BLAST J084546-432458 


263.3005 


-0.2341 


S 


SP&M 


SP 


12 


BLAST J084552-432152 


263.2707 


-0.1886 


S 


SP 


_ 


13 


BLAST J084606-433956 


263.5331 


-0.3425 


s 


M 


_ 


14 


BLAST J0846 12-432337 


263.3322 


-0.1585 


s 


SP&M 


SP 


23 


BLAST J084633-432100 


263.3370 


-0.0827 


P 


SP 


_ 


24 


BLAST J084633-435432 


263.7742 


-0.4309 


P 


SP&M 


SP 


26 


BLAST J084637-432245 


263.3685 


-0.0899 


p 


SP 


_ 


31 


BLAST J084654-43 1626 


263.3177 


0.0149 


s 


SP&M 


SP 


34 


BLAST J0847 18-432804 


263.5147 


-0.0496 


s 


SP&M 


SP 


38 


BLAST J08473 1-435344 


263.8713 


-0.2886 


p 


SP 


SP 


40 


BLAST J084735-432829 


263.5521 


-0.0141 


s 


SP&M 


SP 


41 


BLAST J084736-434332 


263.7488 


-0.1699 


s 


SP&M 


SP 


43 


BLAST J084738-434931 


263.8305 


-0.2277 


p 


SP 


SP 


44 


BLAST J084739-432623 


263.5322 


0.0167 


p 


SP 


_ 


45 


BLAST J084742-434347 


263.7637 


-0.1582 


p 


SP 


SP 


47 


BLAST J084744-435045 


263.8591 


-0.2249 


s 


SP&M 


SP 


50 


BLAST J084749-434810 


263.8349 


-0.1861 


s 


SP&M 


SP 


52 


BLAST J084754-432748 


263.5802 


0.0387 


p 


SP 


_ 


53 


BLAST J084759-433942 


263.7428 


-0.0757 


p 


SP 


_ 


54 


BLAST J084801-435 108 


263.8947 


-0.1907 


p 


SP 


_ 


55 


BLAST J084803-433051 


263.6369 


0.0281 


s 


M 


_ 


56 


BLAST J084805-435415 


263.9430 


-0.2138 


s 


SP&M 


_ 


57 


BLAST J0848 13-423730 


262.9644 


0.6103 


s 


SP&M 


SP 


59 


BLAST J0848 15-4347 14 


263.8713 


-0.1166 


p 


SP 


SP 


63 


BLAST J084822-433152 


263.6856 


0.0609 


s 


SP&M 


SP&M 


65 


BLAST J084823-433858 


263.7799 


-0.0106 


s 


SP 


_ 


71 


BLAST J084834-435455 


264.0059 


-0.1538 


p 


SP 


_ 


72 


BLAST J084834-432430 


263.6126 


0.1657 


s 


SP 


_ 


77 


BLAST J084842-43 1735 


263.5392 


0.2584 


p 


SP 


SP 


79 


BLAST J084844-433733 


263.8007 


0.0525 


p 


SP 


_ 


81 


BLAST J084847-425423 


263.2482 


0.5133 


p 


SP&M 


SP 


82 


BLAST J084848-433225 


263.7415 


0.1155 


p 


SP&M 


SP 


88 


BLAST J0849 10-44 1636 


264.3543 


-0.2984 


p 


SP 




89 


BLAST J0849 12-43 1353 


263.5480 


0.3668 


s 


SP&M 


SP 


90 


BLAST J084912-433618 


263.8379 


0.1312 


p 


SP 


SP 


93 


BLAST J084925-431710 


263.6165 


0.3645 


p 


SP 


SP 


97 


BLAST J084932-441046 


264.3206 


-0.1857 


p 


SP 


SP 


101 


BLAST J084952-433808 


263.9381 


0.2058 


s 


M 




109 


BLAST J085033-433318 


263.9551 


0.3532 


s 


M 





Notes. (a) We follow the source numbering (0 to 140) defined by Olmi et al. (2009). ib) "SP" stands for single-point observation and "M" stands 
for map. 



The parameter 77 m b to convert from antenna temperature to 
main-beam brightness temperatur e has been assum ed to be 0.44 
at 100 GHz and 0.49 at 86 GHz (lLadd et al.ll2005h . The Mopra 
spectrometer (MOPS) was used as a backend instrument in its 
"zoom" mode, which allowed to split the 8.3 GHz instantaneous 
band in up to 16 zoom bands; each sub-band was 137.5 MHz 
wide and had 2x4096 channels. The single-point observations 
were performed in position-switching mode, whereas the spec- 
tral line maps, of size mostly 3x3 arcmin 2 , were obtained using 
the Mopra on-the-fly mapping mode, scanning in both right as- 
cension and declination. 

The most intense lines observed were N2H + (1 -0), HCN(1 - 
0), HNC(1 - 0) and HCO + (l - 0). These lines are all good trac- 
ers of dense gas, and are also often used to detect velocity gradi- 
ents. N2H + is known to be more resistant to freeze-out on grains 
than the car bon-bearing speci es . HNC is particularly prevalent 
in cold gas dHirota et al.ll 19981) . while HCO + often shows infall 
signatures and outflow wings. These strong lines can all be op- 



tically thick and thus two isotopologues, H CO + and H CN, 
were also observed to provide optical depth and line profile in- 
formation. Unfortunately, their intensity was in general too weak 
to allow derivation of column densities throughout the maps and 
other useful physical and kinematical parameters, and thus most 
of our analysis will be based on the four most intense molecular 
tracers. 

2.3. Parkes 

A total of 22 sources where observed with the Parkes tele- 
scope in single-pointing mode (see TableQ]), and only one source 
(BLAST063) was mapped using the NH 3 (1,1) and (2,2) transi- 
tionfl The system temperatures were typically ~ 40 - 50 K and 
we reached a RMS sensitivity in T* units of ~ 10 mK after re- 



1 The signal-to-noise ratio (SNR) in this map was very low and thus 
we will not consider any further this map. 
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binning to a ^ 0.2 kms -1 velocity resolution. The beam FWHM 
was about 1.3 arcmin and the pointing errors were found to be 
< 20arcsec. The Parkes Multibeam correlator, in single-beam 
high resolution mode, was used as a backend. A simultaneous 
band of 64 MHz was observed, with a total of 8192 channels, 
which allowed to observe the NH3(1,1) and (2,2) lines within 
the same bandwidth. 

2.4. BLAST 

BLAST is a 1.8 m Cassegrain telescope, whose under- 
illuminated primary mirror has produced in-flight beams with 
FWHM of 36 ", 42 ", and 60 " at 250, 350, and 500 yum, respec- 
tively. More details on the instrument, calibration methods and 
map-making proced u re can be found i n the BLA ST05 papers 
iPascale et all d2008l) . iPatanchon eTail d2008l) . and iTruch et all 
(120081) . 

The Vela-D r egion, shown in Figure Q] and defined by 
lOlmi et all d2009l) as the area contained within 262?80 < I < 
264!60 and -1°15 < b < 1°10, was part of a larger map 
of the Galactic Plane obtained by BLAST toward the VMR 
dNetterfield et al.ll2009l) . In Vela-D lOlmi etal] d2009l) found a to- 
tal of 141 sources, both starless and proto-stellar, and in the next 
sections we will follow the source numbering (0 to 140) defined 
by these authors. 

3. Results 

In this section we will discuss the derivation of various physical 
parameters such as column density, mass and kinetic tempera- 
ture, that are fundamental for the analysis of the sources. We 
will also analyse the kinematics of the sources, determining ve- 
locity gradients and line asymmetries. The estimate of the virial 
masses will also allow us to determine which sources are cur- 
rently gravitationally bound and which are not. 

3. 1 . Morphological characteristics of spectral line maps 

The spectral line data and maps were reduced and analyzed us- 
ing the standard IRAM package CLASS, and the XS package of 
the Onsala Space Observatory. The single-point spectra obtained 
with the Mopra telescope towards the nominal position of the 
BLAST cores are shown in Figures |B~T1 to IB ,7l of the Appendix. 
The peak intensities of the detected lines, or the spectrum RMS 
in case of no detection, are listed in Table [3] Tables ICTI to IC. 41 
in the Appendix list several physical and kinematical parameters 
of the sources. 

The multi-line integrated intensity maps obtained at Mopra, 
with the most significant SNR, are shown in Figure |A~T1 to I A.4| 
of the Appendix, where the white contours represent the BLAST 
emission at 250 fim. One can see that in general the line emission 
follows the dust continuum emission, even when the emission is 
more diffuse and filamentary as observed, for example, in source 
BLAST063. But we also note that in several sources (for exam- 
ple, BLAST055, BLAST056, BLAST08 1 and BLAST101) the 
dust continuum and some of the line emission peak at different 
positions. 

One can immediately note that the three proto-stellar cores 
that were mapped (BLAST024, BLAST081 and BLAST082) 
have a much more regular shape than the starless ones. We also 
note that the maps of sources BLAST024 and BLAST081 are 
partly affected by artefacts, likely caused by bad weather during 
the observations at Mopra, which show up as sharp drops of the 



Table 4. NH 3 (1,1) line parameters. 



Source # 




NH 3 (1,1) 








Vlsr 


FWHM 


r 




[K] 


[kms" 1 ] 


[kms" 1 ] 




3 


0.15 


5.3 


1.9 


1.40 


9 


0.04 


10.7 


1.1 


0.66 


14 


0.13 


3.0 


0.7 


2.63 


24 


0.38 


3.5 


2.6 


0.51 


31 


0.06 


10.7 


1.3 


2.77 


40 


0.07 


1.3 


0.9 


0.22 


41 


0.05 


11.1 


0.7 


0.10 


47 


0.04 


5.5 


1.5 


0.77 


57 


0.07 


14.0 


1.5 


1.73 


63 


0.11 


12.4 


1.7 


0.94 


77 


0.05 


2.0 


0.8 


0.10 


81 


0.06 


12.3 


3.1 


0.10 


82 


0.11 


11.7 


2.0 


0.22 


90 


0.10 


11.4 


1.7 


0.10 


93 


0.13 


2.4 


1.2 


0.70 


97 


0.31 


9.4 


1.6 


1.28 



Notes. Parameters derived from method "nh3" of the CLASS program. 
Sources where the hyperfine structure of the NH 3 (1,1) could not be fit- 
ted are not listed here. 

Table 5. NH3(2,2) line parameters. 



Source # 




NH 3 (2,2) 








v Isr 


FWHM 




[Kkms- 1 ] 


[kms" 1 ] 


[kms" 1 ] 


3 


0.06 


5.3 


1.1 


24 


0.54 


3.6 


3.1 


47 


0.03 


5.3 


2.9 


63 


0.03 


10.3 


1.8 


81 


0.08 


12.2 


2.7 


82 


0.07 


11.7 


1.9 


90 


0.06 


11.4 


2.5 


93 


0.07 


2.4 


1.8 


97 


0.14 


9.3 


2.1 



Notes. Parameters derived from a standard Gaussian fit to the main 
component of the NH 3 (2,2) hyperfine structure. 



integrated intensity, particularly in the maps of the HCO + (l - 0) 
line. We can also see that in source BLAST082 two separate, 
smaller cores have been identified. In terms of the molecu- 
lar tracers observed, we note that there are clear differences 
in the spatial distribution of emission from different molecular 
species. These morphological differences are more evident in 
the starless cores (but also, for example, in the proto-stellar core 
BLAST08 1) and suggest both physical and chemical differences 
among the observed sources. 

We may summarize the wide variety of spatial morpholo- 
gies observed as follows: very early and cold starless cores ap- 
pear to have an irregular shape (e.g., BLAST009) in most or 
all molecular tracers mapped at Mopra. Some warmer cores 
(e.g., BLAST031) and cores at the transition phase from star- 
less to proto-stellar (BLAST063, see Section |43T > appear to be 
more compact. Finally, proto-stellar cores all show a more reg- 
ular shape and narrow radial intensity profiles (see Section [3~6l >. 
Throughout the rest of the paper we will thus discuss how, be- 
sides to morphological differences, these cores also show physi- 
cal and chemical variations, which may be associated to different 
evolutionary phases. 
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Table 3. Peak intensity of main detected lines or the rms Level of the spectrum. 



Source # 


N 2 H + 


HCN 


H 1J CN 


HNC 


HC 3 N fl 


HCO + 


H 1J CO + 


3 


0.17 


0.81 


<0.03 


0.58 


- 


0.70 


0.13 


9 


0.09 


0.06 


<0.03 


0.20 


- 


0.22 


0.06 


12 


<0.02 


0.05 


<0.02 


<0.03 


- 


0.07 


<0.02 


23 


<0.02 


0.10 


<0.03 


0.25 


- 


0.38 


<0.03 


24 


1.33 


2.97 


0.23 


2.23 


- 


3.20 


0.32 


26 


<0.02 


0.07 


<0.02 


0.08 


- 


0.14 


<0.02 


26'' 


- 


- 


- 


0.10 


- 


0.13 


- 


31 


<0.02 


0.11 


<0.02 


0.15 


- 


0.34 


<0.02 


34 


<0.02 


0.08 


<0.02 


0.07 


- 


0.23 


<0.02 


34 c 


- 


- 


- 


- 


- 


0.15 


- 


38 


0.07 


0.11 


<0.02 


0.15 


- 


0.19 


<0.02 


40 


0.20 


0.31 


<0.02 


0.45 


- 


0.73 


0.28 


41 


<0.02 


0.27 


<0.03 


0.20 


- 


0.41 


<0.02 


43 


<0.05 


0.11 


<0.05 


0.26 


- 


0.42 


<0.05 


44 


0.21 


0.35 


<0.04 


0.57 


- 


0.88 


0.13 


45 


0.13 


0.24 


<0.02 


0.24 


- 


0.52 


<0.02 


47 


<0.02 


0.19 


<0.02 


0.12 


- 


0.14 


<0.02 


47'' 


- 


- 


- 


- 


- 


0.13 


- 


50 


<0.02 


0.10 


<0.02 


0.07 


- 


0.16 


<0.02 


52 


<0.02 


0.11 


<0.02 


0.10 


- 


0.18 


<0.02 


53 


0.07 


0.20 


<0.02 


0.28 


- 


0.19 


0.07 


54 


<0.02 


0.14 


<0.02 


0.12 


- 


0.12 


<0.02 


56 


<0.02 


0.14 


<0.03 


0.08 


- 


0.16 


<0.02 


57 


0.12 


0.13 


<0.03 


0.21 


- 


0.33 


0.13 


59 


0.12 


0.18 


<0.03 


0.22 


- 


0.23 


0.08 


63 


0.11 


0.44 


<0.02 


0.55 


- 


0.76 


0.09 


65 


<0.02 


0.12 


<0.02 


0.09 


- 


0.15 


<0.02 


71 


<0.04 


0.11 


<0.04 


<0.05 


- 


<0.08 


<0.04 


72 


<0.02 


0.15 


<0.02 


0.14 


- 


0.39 


<0.02 




- 


- 


- 


- 


- 


0.19 


- 


77 


0.06 


0.26 


<0.02 


0.23 


- 


0.45 


<0.03 


79 


<0.05 


0.17 


<0.05 


0.23 




0.22 


<0.05 


81 


1.49 


4.32 


0.39 


2.42 




4.31 


0.25 


82 


0.50 


0.65 


<0.03 


0.74 




1.30 


0.23 


88 


<0.02 


0.26 


<0.03 


0.18 




<0.02 


<0.03 


89 


<0.02 


0.16 


<0.02 


0.09 




0.22 


<0.03 


90 


0.50 


1.32 


0.12 


1.30 


0.16 


2.19 


0.38 


93 


0.50 


0.75 


0.11 


1.03 


0.26 


1.45 


0.44 


97 


0.19 


0.54 


<0.03 


0.62 




0.79 


0.12 



Notes. Only sources with a single-point observation are listed (see Table[T]l. All values are in units of T£ [K]. 

(a) The typical noise RMS in the spectra where HC3N was not detected is a 0.03 K. ™ This source has two separate velocity components in 
HNC and HCO + and only these lines are listed for this velocity component. <rt These sources have two separate velocity components in HCO + 
and only this line is listed for this velocity component. 



3.2. Temperature and optical depth from hyperfine structure 
fits 

The excitation temperature, T ex , and the line optical depth, t, 
which are required to estimate the molecular column densi- 
ties (see Sect. 13.5. U could only be determined in those sources 
where molecules with a hyperfine structure (hfs, hereafter) were 
detected, i.e., N 2 H + , HCN and H 13 CN, and/or when both the 
HCO + (l - 0) and H 13 CO + (l - 0) lines were detected. 

In the molecules with a hfs, we used method "hfs" of 
the CLASS program to determine both T ex and t. In those 
cases where both HCO + and H 13 CO + were detected, the opti- 
cal depth was estimated from their line intensity ratio, assuming 
the same excitation temperature for the two isotopic molecular 
species and using the abundance ratio [HCO + ]/[H 13 CO + ]- 40 
dZinchenko et alj|2009l) . 

For the purpose of estimating the column density later, in 
those sources where the hfs method could not be used for one 
or more molecules (for example because of low SNR), or if only 



HCO + was detected but not H 13 CO + , the values of both T ex and t 
for a given molecular species were assumed to be the averages of 
the parameters estimated from the detected molecular transitions 
in the same source. 

3.3. Temperature and optical depth from NH 3 

The single-point spectra of the NH3(1,1) inversion line, ob- 
tained with the Parkes telescope, are shown in Figure lB~8l of the 
Appendix. The (2,2) line has been detected towards nine sources 
(mostly proto-stellar). We attempted to map the NH3(1,1) tran- 
sition toward BLAST063, but we did not detect the line due to 
insufficient SNR in the individual scanning positions. 

Method "nh3" of the CLASS program, similar to method 
"hfs" described above, was used to fit the NH3(1,1) line and re- 
sulted in the parameters T% t, t and FWHM listed in Table [4] 
Due to the fact that the hfs components of the NH3(2,2) line 
were not detected (with the exception of source BLAST024), 
the line parameters were derived from a standard Gaussian fit to 
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Table 6. Results of velocity gradient fitting. 



Source # 




HNC(l-O) 






HCO + (l -0) 






V 


dV/dr 


6/ 


V 


dV/dr 


e v a 




[kms 4 ] 


[km s per 1 ] 


[deg] 


[km s~'] 


[kms -1 per 1 ] 


[deg] 


9 


- 


- 


- 


12.2 


0.8 


88.5 


13 


1.0 


2.3 


21.1 


0.9 


1.5 


18.1 


24 


5.5 


1.9 


60.7 


5.3 


1.7 


57.7 


34 








2.6 


1.6 


224.1 


40 


2.5 


0.2 


16.7 


2.6 


1.2 


71.5 


41 








12.1 


1.2 


289.9 


47 


6.6 


2.7 


240.1 








63 


13.0 


1.6 


55.2 


13.2 


1.3 


86.8 


81 


13.5 


1.4 


256.9 


13.3 


1.9 


312.6 


82 


13.6 


2.1 


15.7 


13.0 


1.5 


334.5 



Notes. The angle V is measured positive from the axis of positive longitude offsets toward the North. 
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Fig. 2. Kinetic temperature, 7\, as derived from the NH3(1,1) 
and (2,2) inversion transitions, vs. the BLAST-derived dust 
temperature, Sources are identified through their ID (see 
Table [TJ|, and the vertical arrows indicate that the correspond- 
ing 7\ is an upper limit. The dashed line indicates the 77k = T& 
locus. 

the main component of the NH3(2,2) hyperfine structure and are 
listed in Table [5] 

We then determined the kinetic temperature, 7k, and NH3 
column density using both NITi(l,l) and (2,2) transitions, in 
those sources where the (2,2) line was indeed detected. In these 
cases we were able to determine the rotational temperature, 
T n (and thus 7\), and the col umn density us i ng the method 
of lUngerechts et all d!986l) and iBachiller et ail dl987b . In this 
method the required data were: (i) the product r(77 ex - 77b g ) for 
the (1,1) line, where 77 ex and T\, g = 2.725 K are the excitation 
and background temperatures, respectively, and r is the optical 
depth; ( ii) the linewidth A V( 1,1); and f Hi) the (2,2) integrated in- 
tensity. Once Tu has been estimated, the kinetic temp erature can 
be de termined using the analytical expression of Tafalla et al. 
(2004}: 



l-^ln[l + l.lexp(-i|)] 

which is an empirical expression that fits the T^-T^ relation ob- 
tained using a radiative transfer model. For those sources where 



the NH3(2,2) transition was not detected, we could only deter- 
mine an upper limit to 77k (shown in Figure 0. However, in or- 
der to obtain an estimate of the column density, and thus of the 
mass, for these sources we adopted the BLAST-derived temper- 
ature, T&, as an approximation for 7\, which turns out to be a 
reasonably good assumption, as discussed below. 

For those sources with an independent estimate of 7k we 
compared the resulting kinetic temperat ures with the BL AST- 
derived dust temperatures, as obtained by Ol mi et all (2009), and 
we show the results in Figure[2] This comparison between 7\ and 
Td, when applied separately to starless and proto-stellar cores, is 
affected by the few detections of NH3(2,2) toward starless cores, 
and by the relatively high spread in the values of the Tk/Td ratio. 
In fact, for the two starless cores (BLAST047 and BLAST063; 
see also Section l4~5l l where this ratio could be measured we ob- 
tain a mean value (Tk/T&) = 1.4, whereas for the proto-stellar 
cores we get a median value T^/T^ = 1.2 + 0.2. 

When all cores are considered together, we get a median 
value 7k/77d = 1.2 + 0.2, i.e., the median value is dominated by 
the proto-stellar cores. A tentative explanation for the fact that 
7k tends to be slightly higher than 7d could be that the NH3 
observations are sampling regions that are somewhat warmer 
compared to the BLAST-derived T& measurements. In fact, the 
BLAST observations are likely to average the temperature on 
much larger volumes of dust and, in addition, in the less dense 
regions there could be a systematic difference between 7^ and 
7V In the literature both cases of systematic departures between 
dust and gas temperature, i.e. 7\ < or 7\ > 7^, can be foun d 
(e.g.. iKruegel & Walmslevlll984 iPirogov & Zinchenko|[T99l . 
The dust and gas temperature may be affected by various effects, 
and a more detailed analysis of these effects in our sources is 
beyond the scopes of the present work. 

3.4. Kinematics 

3.4.1 . Velocity gradients 

We investigated the presence of systematic velocity gradients in 
those sources which could be mapped. In particular, the velocity 
gradient of a gas clump can be determined by using all or most 
of the data in a map at once, by least-squares fitting maps of line- 
center velocity for t he direction and mag nitude of the best-fit ve- 
locity gradient (see lGoodman et al.ll993l and references therein). 

A cloud undergoing solid-body rotation would exhibit a lin- 
ear gradient, dV/dr, across the face of a map, perpendicular to 
the rotation axis. Thus, we have applied the procedure developed 
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Offset [rod] Offset [red] 

Fig. 3. In the top panels we show the maps of the Vi sl for the HNC( 1-0) (left) and HCO + ( 1 - 0) (right) lines toward the proto-stellar 
core BLAST081, as obtained from Gaussian fits to the spectra at each map position. The arrows show the direction of the velocity 
gradient, with their length proportional to the magnitudes listed in Table|6] The bottom panels show the position-velocity plots, with 
the horizontal axis representing the angular offset along a direction parallel to the velocity gradient shown by the arrows in the top 
panels. Note the different direction of the velocity gradient in the two molecular tracers. 



bv lGoodman et alJ dl993h and lOlmi & Testil d2002h to determine 
the magnitude and direction of the velocity gradient in our spec- 
tral line maps. Our best results are listed in Table [6] and are re- 
stricted to the HCO + (l - 0) and HNC(1 - 0) lines. We note that 
for most sources with measured velocity gradients in both lines, 
the two separately estimated values of dV/dr and 6 V are con- 
sistent, within a 50 % difference. However, the velocity gradient 
measured in the HNC map of core BLAST040 is quite small and 
thus its direction can be affected by large uncertainties. Thus the 
largest discrepancies between the two tracers are observed in the 
two proto-stellar cores BLAST081 and BLAST082, and also in 
the starless core BLAST013. 

The position-velocity plots for the proto-stellar core 
BLAST08 1 are shown in Figure[3] The top panels show the maps 
of V lsr for the HNC(1 - 0) and HCO + (l - 0) lines, as obtained 
from Gaussian fits to the spectra at each observed position. The 
arrows represent the direction and magnitude of the velocity gra- 
dient as obtained from the procedure described above. We note 
the significant difference between the directions of the velocity 
gradients as measured in the HNC and HCO + maps. It is not 
clear whether this discrepancy is real or may be caused by the 
presence of some artefacts in the HCO + map of BLAST081, as 
noted in Section [3~T1 On the other hand, source BLAST082 has 
no visible artefacts, and it is thus possible that in both proto- 
stellar cores the different velocity gradients in the HNC and 
HCO + maps are real and are tracing different systematic mo- 
tions. One of these velocity gradients could be associated with a 



molecular outflow, as suggested by the presence of line- wings in 
the HCO + (l - 0) spectra (see Table IcA. 

3.4.2. Line asymmetries 

The HCO + (l - 0) line is generally the most intense transition 
observed with the Mopra telescope in each source, and it does 
not have a hyperfine structure. Therefore, it represents the ideal 
candidate to identify and analyze asymmetries of the line pro- 
file. In fact, the HCO + (l - 0) line shows a non-Gaussian profile 
toward several sources (see Table lC.4l i. and in some cases unam- 
biguous wing emission is detected, indicating the likely presence 
of a molecular outflow. By inspecting Table IC.41 we also note a 
moderate shift, up to ~ 1 kms -1 , between the peak positions of 
the generally optically thick HCO + (l - 0) line and optically thin 
H 13 CO + (l - 0) transition. We can give a quantitative estimate of 
this asymmetry in order to get some indirect information about 
dynamical processes in the cores. 

A widely used method to extract line asymmetries is based 
on the comparison of optically thin and optically thick line 
velocities, i. e., determining the p arameter 5V = (V t hick - 
Vthin)/AVthin dMardones et al][l997D . In this way we can divide 
our cores in "blue shifted" cores, with 6V < 0, and "red shifted" 
cores with 6V > 0. The blue shifted values of 6V could be caused 
by infall motions and the red excess by expanding motions or 
outflow. 

We identified 13 sources where the HCO + (l - 0) and 
H 13 CO + (l - 0) linewidths could be measured and in Figure [4] 
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we show the distribution of their 6V parameter, with V tn i c k and 
Vthin being determined from the HCO + (l -0) and H 13 CO + (l -0) 
lines, respectively. We also estimated the typical error on 5V and 
then used the same +5cr$y threshold as M ardones et al.l {l991) 
to exclude those sources where the Vi sr differences were domi- 
nated by measurement errors. We estimated 5cr sv - 0.4 - 0.6, 
thus in Figure [4] a typical value 5cr 5V =2 0.5 is used. Only three 
cores show a significant blue shift (BLAST044, BLAST090 and 
BLAST093), and two a red shift (BLAST003 and BLAST053). 
From the red excess candidates we have to exclude BLAST053 
which has a weak and noisy H 13 CO + (l - 0) spectrum. 

3.5. Derivation of masses 

From our spectral line data we can derive masses from the gas 
column density, M c d, as well as virial masses, M vu -, which we 
can then compare with the BLAST-derived masses, Mbiast, esti- 
mated from the dust continuum. Clearly, while the virial masses 
refer to the total gas mass, with M cd we can only determine 
the mass of a given molecular species. The total gas mass can 
only be determined by assuming specific molecular abundances; 
or, alternatively, we can give a coarse estimate of the molecular 
abundance by taking the ratio M c< j/Mbiast (see Section |4~T| |. The 
isotopologues observed by us, H 13 CO + and H 13 CN, are gener- 
ally too weak to allow derivation of column densities throughout 
the maps, and thus we use only the main molecular species, at- 
tempting to correct for first-order optical depth effects. 

3.5.1 . Mass derived from column density 

The calculation of masses from the column density assumes 
LTE and we determine the molecular gas mass integrating the 
molecule column density over the extent of the source. Then we 
can write: 



M cd = d 2 m moi N mol dO 



(2) 



where J N mo \ dQ is the molecule column density integrated over 
the region enclosed by the chosen contour level, m mo \ is the mass 
of the specific molecule being considered, and d is the distance 
to the source. The column density N mo \ corresponding to a J — » 
7-1 rotation transition can be calculated as: 



AUi [cm 2 ] = 
-— f 



4.0 x 10 12 
7V[D] S[K] 



Z exp — 

\ * ex 



r*dv [Kkms- 1 ] 



(3) 



where B denotes the rotational constant, £j is the upper state 
energy, // is the dipole moment (in Debye) and we used the es- 
cape probability t/[1 - exp(-r)] to account for first-order op- 
tical depth effects. The derivation of both T ex and r has been 
discussed in Section [3~!2l For kT ex » hB the partition function, 
Z, of a linear molecule is given by: 



~~hB~ 



(4) 



where h and k are the Planck and Boltzmann constants, respec- 
tively. 

Eq. (f2j) is actually implemented by writing: 



M cd = d 2 m mo[ Afi pix 2^ N mo i(Xi,yi) 



(5) 



6 h 



c 

m 5 



CD 



4 h 
3 r 
2 r 
1 \ 

L 



1 .0 



-0.5 



0.0 

dV 



0.5 



.0 



Fig. 4. Histogram of 6V = (Vthick- Vthin)/ A Vthin, where V th i c k and 
Vhin were determined from the HCO + (l - 0) and H 13 CO + (l - 0) 
lines, respectively. The vertical dashed lines indicate the ±5o- sv 
level. 



where N mo \(xi,yi) represents the column density in a single pixel 
(x,-,y,) of the map, with n p ; x representing the total number of 
pixels, and AQ p ; x represents the solid angle covered by a sin- 
gle pixel. The map pixels selected are those that have an in- 
tegrated intensity / > cr map , and also lie within an 1/e radius 
from the center of a 2D Gaussian fit to the core, for each specific 
molecule. Assuming that a line is detected if at least two adjacent 
velocity channels lie above the 3<T rms level of the spectrum, then 
we calculate the RMS of the integrated intensity in the map sim- 
ply as cr map = 2(3cr rms )Av, where Av is the channel width. We 
thus select only those pixels that belong to the actual "core" of 
the source, and remove pixels that come from extended diffuse 
envelopes. The resulting masses are listed in Table|7] 



3.5.2. Virial mass 

We then estimate the virial mass of the cores as suming they are 
simpl e spherical systems with uniform density dMacLaren et al.l 



M vir [M ] = 2107? dec [pc] (AV av [km s" 1 ]) 2 



(6) 



where is the deconvolved source radius and AV av represents 
the line FWHM of the molecule of mean mass, calculated as the 
sum of the thermal and turbulent components: 



A V 2 , = AV 2 ol + kT8 In 2 j — — 

\m av m mo i 



(7) 



where AV mo i is the line FWHM of the molecular transition 
being considered and m av = 2.3 amu is the mean molecular 
weight (with respect to the total number of particles), assum- 
ing a mass frac t ion fo r He of 25%. We note that according to 
iMacLaren et al.l d!988l) . Eq. (O may lead to a mass overestimate 
if the density distribution is not uniform. For example, in the 
case of a sphere with a density distribution p oc r~ 2 the numeri- 
cal factor in Eq. (0 should be replaced by 126. 

As the decon volved sour c e radi us we used the BLAST- 
derived values of Olmietal ] (I2009h . We decided to use the 
BLAST sizes, to estimate the virial masses, instead of the 
Mopra-derived sizes for two main reasons: (i) the linewidths 
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Table 7. Mass estimates. 



Source # 


















M Ua8 , 


M vil fl 


NiH + 


HNC 


TT7T7TT 

HCO + 


HCN 


NH 3 C 




[M Q ] 


[M ] 


[xl0- 9 M G ] 


[xlO- 5 M G ] 


[xlO- 9 M Q ] 


[xl0- 9 M o ] 


[xl0- 7 M o ] 


3 


23.2 


68.1 


— 


— 


— 


— 


— 


9 


3.1 


15.9 


— 


2.6 


1.6 


0.3 


7.9 


13 


6.4 


— 


— 


3.7 


7.6 


3.5 


— 


14 


7.7 


8.9 


— 


— 


— 


— 


— 


24 


98.0 


87.3 


5.4 


25.7 


54.1 


17.8 


5.8 


31 


4.0 


26.4 


— 


0.02 


2.3 


0.6 


31.8 


34 


4.6 


— 


— 


2.1 


1.9 


— 


— 


38 


2.2 


32.4 


— 


— 


— 


— 


— 


40 


6.3 


11.3 


0.13 


4.9 


13.9 


2.5 


2.1 


41 


2.8 


10.8 


— 


— 


7.4 


— 


0.8 


44 


5.8 


19.1 


— 


— 


— 


— 


— 


45 


2.1 


17.6 


— 


— 


— 


— 


— 


47 


2.6 


44.7 


— 


2.2 


2.7 


5.7 


4.6 


50 


2.3 


— 


— 


— 


6.0 


9.9 


— 


53 


11.6 


25.9 


— 


— 


— 


— 


— 


56 


1.8 


— 


— 


0.1 


1.5 


0.6 


— 


57 


14.0 


27.0 


— 


— 


8.7 


— 


— 


59 


1.8 


17.5 


— 


— 


— 


— 


— 


63 


13.4 


64.6 


1.5 


4.7 


10.3 


8.3 


9.0 


77 


3.7 


12.3 












81 


6y.y 


oo o 
88.8 


2.2 


1 Q A 

1.3.4 


/in a 
40.4 


1 A 1 

14.2 


U.O 


82 


16.5 


42.2 


2.4 


8.1 


29.8 


4.2 


1.4 


89 


2.4 








6.7 






90 


22.2 


43.9 












93 


15.2 


17.7 












97 


24.4 


54.5 












101 


2.1 






1.0 


3.5 


0.9 




109 


1.0 








0.9 







Notes. <a> M v j, is determined assuming a uniform density spherical source. <b> is determined only when a spectral line map was available, 
with the exception of NH 3 . u) M[NH 3 ] is listed only for those sources where an estimate of the rotational temperature, T i2 , was possible, and 
where an estimate of the source diameter could be obtained from at least one of the molecular tracers. 

are mostly measured from the single -point spectra of NH3, and 
where NH3 was not observed or not detected we used the ob- 
servations of N2H + , for which only a few maps were available; 
( ii) the BLAST maps are not affected by spatial variations that 
may instead affect specific molecules and thus likely give a bet- 
ter representation of the mass distribution in each source. 

As previously stated, the line FWHM is calculated from the 
NH?(1,1) single-point spectra toward each source, if available, 
or from the N2H + (1 -0) line otherwise, since these molecules are 
less likely to be affected by depletion, they trace the denser gas 
and are mostly optically thin. The linewidths determined by the 
fit procedure are artificially broadened by the velocity resolution 
of the observations, and thus we subtracted in quadrature the 
resolution width, AV, es , from the observed line FWHM, AV t, s , 

such that AV raol = ^AV o 2 bs - AV* S . 

The resulting virial masses are listed in Table [7] and a plot 
of the total core mass, Mbiast, as derived from BLAST, vs. the 
virial mass is shown in Figure [5] We note that nearly all of 
the starless cores lie below the self-gravitating line, indicating 
that they are unlikely to be gravitationally bound, whereas more 
than half of the proto-stellar cores lie near or above that line, 
confirming that they are gravitationally bound (about 30% of 
all sources have M core /M v i r > 0.5). However, as we mentioned 
above, for a non-uniform density distribution (see Section [3~6l l. 
e.g. of type p oc r~ 2 ', all virial masses should be multiplied by 
a factor 0.6. In this case, most of the proto-stellar cores would 
move above the self-gravitating line and also most of the starless 



100 r 




10 



100 



M vir [M ] 



Fig. 5. Total core mass, M core = Mbi as t, vs. the virial mass, M v - lr , 
calculated using the velocity linewidths of the NH3(1,1) and 
N2H + (1 - 0) transitions (see text). Yellow, filled squares repre- 
sent proto-stellar cores and red, filled triangles represent starless 
cores. The solid line indicates the minimum M core = 0.5M v i r for 
which the cores should be self-gravitating. Individual values of 
M v j r are listed in Table|7] 
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BLAST063 BLAST063 BLAST063 BLAST063 




02468 02468 02468 02468 

Ring Radius [pixels] Ring Radius [pixels] Ring Radius [pixels] Ring Radius [pixels] 



Fig. 6. Examples of normalized radial profiles for sources BLAST063 (starless, top panels) and BLAST082 (proto-stellar, bottom 
panels). The points and error bars represent the ring-averaged integrated intensity in the molecular lines (from left to right) N2H + ( 1 - 
0), HNC(1 - 0), HCO + (l - 0) and HCN(1 - 0). The black solid line represents the fit obtained using Eq. ®. The red solid line 
shows the radial profile of the averaged BLAST intensity at 250/im, and the green solid line represents the 3-mm beam profile of 
the Mopra telescope. Each pixel corresponds to 15 arcsec. 



Table 8. Parameters for fit to the radial profile of column density. 



Source # 


N 2 H + 




HNC 




HCO + 




HCN 




Rflat 

[arcsec] 


1 


Rflat 

[arcsec] 




Rflat 

[arcsec] 




Rflat 

[arcsec] 


1 


9 






67.5 


1.2 


15.0 


0.3 


34.5 


0.9 


13 






69.0 


3.4 


16.5 


0.6 


19.5 


0.5 


24 


28.5 


1.1 


48.0 


1.5 


15.0 


0.8 


22.5 


1.1 


31 






37.5 


3.3 


27.0 


1.8 


15.0 


0.4 


34 






57.0 


1.4 


16.5 


0.7 






40 


27.0 


2.7 


45.0 


1.0 


16.5 


0.3 


16.5 


0.2 


41 










33.0 


0.8 






47 






30.0 


0.8 


16.5 


0.4 


15.0 


0.4 


50 










15.0 


0.3 


18.0 


0.6 


56 






58.5 


2.3 


39.0 


1.9 


16.5 


0.5 


57 










20.0 


0.6 






63 


70.0 


2.0 


30.0 


0.9 


18.0 


0.4 


24.0 


0.5 


81 


24.0 


2.6 


21.0 


1.3 


15.0 


0.7 


18.0 


1.1 


82A" 


31.5 


2.4 


22.5 


0.7 


16.5 


0.5 


21.0 


0.5 


82B 


31.5 


2.4 


16.5 


0.6 


16.5 


0.5 


16.5 


0.5 


89 










36.0 


1.1 






101 






22.5 


0.7 


27.0 


1.2 


33.0 


1.6 


109 










36.0 


2.8 







Notes. (fl) The parameters for the two spatial components in BLAST082 are listed. 



cores would be located near that line (about 60% of all sources 
have M core /M v j r > 0.5 in this case). 



3.6. Radial profiles of column density 

In order to identify significant differences in the distribution of 
both molecular gas and dust continuum, we have derived radial 
profiles of the column density from both the BLAST maps at 
250 fim and from the Mopra maps of the integrated intensity of 
the different molecular transitions. While the BLAST 500 /im 
measurements would be somewhat less affected by optical depth 
effects, we have chosen to compare our molecular data to the 
BLAST 250 //m measurements since they are closely matched 
to the Mopra beam. 

Previous work has shown that single-power-law density, or 
column density, profiles do not fit the emission from dense cores 



and that a central flattening is always needed to reproduce the 
data (e.g., Taf alla et al.l l2002 and references therein). Among the 
standard analytic profiles that combine the power-law behavior 
of column density for large radii r and a central flattening at 
small r, is the following function: 



N{r) = /Y fla 



^fla 



flat 



. r 2)l/2 



(8) 



where the column density is approximately uniform, with N ~ 
Nflat> f°i" r ^ ^flat, and it falls off as r~ n for r » R^. 

In the fitting procedure we first circularly average the maps 
of the integrated intensity of selected molecular transitions 
around the peak, and then use a% 2 routine to fit the column den- 
sity profile, obtained from the integrated intensity with Eq. (0, 
using the model given by Eq. ([8]), after convolution with a 2D 
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Table 9. Average properties of radial profiles. 



Spectral line 


(Ram) 
[arcsec] 


(v) 




(Xmol) 




N 2 H + (1 -0) 


35.4 ± 17.2 


2.24 ± 0.65 


(7.3 


+ 5.4) x 10" 


ii 


HNC(l-O) 


40.4 ± 18.3 


1.39 ± 1.00 


(4.5 


± 2.9) x 10- 


10 


HCO + (l -0) 


21.9 + 8.6 


0.87 ± 0.67 


(1.3 


± 0.8) x 10" 


-9 


HCN(1 -0) 


20.8 + 6.4 


0.68 ± 0.39 


(0.8 


± 1.2) x 10" 


-9 


NH 3 a 






(1.7 


±3.1) x lO- 


-7 



Notes. <a) Mean value estimated from NH 3 abundances listed in 
Table [H 



4. Discussion 

In this section we will describe how the results obtained in 
Section [3] allow us to determine the relative molecular abun- 
dances, and in particular how these abundances vary as a func- 
tion of the radial distance in both starless and proto-stellar cores. 
We will also discuss how the relative molecular abundances can 
be used to roughly estimate the chemical evolutionary status of 
the sources. In the last part of the section we will analyze the 
degree of correlations among various physical and chemical pa- 
rameters, which will also allow to highlight further differences 
between starless and proto-stellar cores. 



38" Gaussian. The procedure is then repeated for the BLAST250 
maps. The column density profiles obtained from the spectral 
line maps and the dust continuum maps are then compared to 
find any possible evidence of chemical effects such as, for ex- 
ample, freezing-out on grains. 

The results of the radial profile fits to individual sources are 
listed in Table [8] while in Table [9] we list the average values of 
the parameters rj and /?fl at , i.e. (rj) and (Rfl at ), obtained by se- 
lecting only fits to the radial profiles that converged. We can 
note several trends: first, apart from the large scatter observed 
in the case of the N2H + and HNC molecules, the (Rfl at ) val- 
ues are quite consistent within the errors. Second, the value of 
(rj) for the N2H + molecule is significantly higher compared to 
the other molecules although, once again, the HNC molecule 
shows a larger scatter. We tentatively note a trend of decreas- 
ing (rj) values from HNC to HCO + and then HCN, but given the 
errors this cannot be considered significant. We also note that 
the large uncertainty in the (Raat) value for N2H + is due to the 
BLAST063 source. If this source is not included, then we would 
have (Rflat[N2H + ]) = 28.5 + 3.2, quite close to the values esti- 
mated for HCN and HCO + . 

As a final trend, when comparing the radial profiles obtained 
by the BLAST dust continuum with those of the molecular trac- 
ers, Figure [6] shows that, at least for the two cores shown in the 
figure, the agreement between the radial profiles of the molec- 
ular tracers and of the dust continuum is better for the starless 
core (BLAST063). In order to get a quantitative comparison, 
we have applied the two-sample Kolmogorov-Smirnov test (K-S 
test) to the BLAST- and Mopra-derived radial profiles. The K- 
S test evaluates the maximum deviation between the cumulative 
distribution functions of the two data sets and also the associ- 
ated probability, p^ s (0 < p^ s < 1), that two arrays of data values 
are drawn from the same distribution. Small values of p^ show 
that the cumulative distribution function of one data set is sig- 
nificantly different from the other. For the two sources shown in 
Figure [6] we get an average value of p^ = 0.66 for BLAST063 
and = 0.24 for BLAST082, confirming that in BLAST063 
the dust and molecular gas radial profiles are more similar than 
in BLAST082. 

However, it is clear comparing the radial profiles and also 
the integrated intensity maps, that the starless cores have a wider 
range of morphologies, with sources showing a good agreement 
between the dust continuum and the molecular radial profiles (al- 
though not necessarily in all molecular tracers, e.g., in sources 
BLAST013 and BLAST031), and sources where this agreement 
is less evident (e.g., in BLAST056 and BLAST063). These dif- 
ferences suggest that the starless sources in our sample are cur- 
rently in different evolutionary phases. 



4.1. Determining molecular abundances 

As we mentioned earlier, we have estimated the molecular abun- 
dances by comparing the molecular gas mass, M c( j, with the total 
gas mass, Mbi ast , estimated by the BLAST measurements of the 
dust continuum, i.e., X mo \ = M c d/Mbi ast . We prefer this method to 
the alternative technique of obtaining X mo \ by means of a pixel- 
by-pixel comparison of the column densities in the Mopra and 
BLAST maps, which could be affected by differences in the spa- 
tial distribution of the spectral line and dust continuum emission. 
The estimated abundances are shown in Table \10\ and their av- 
erage values are listed in Table [9] The large uncertainties in the 
estimated abundances are likely the consequence of uncertainties 
in the derivation of the BLAST masses, of the different spatial 
distribution of the molecular gas and the dust, and possibly also 
of different chemical evolutionary phases (see the discussion in 
Sections l4~2landl4~3l. 

The average value of X(N 2 H + ) = (7.3 + 5.4) 1Q-" listed 
in Tab le [9] is consistent with that found by, e.g.. | Blake et al.l 
(I1995I) in NGC1333 and also bv lZinchenko et al.l(l2009l) in SI 87 
and W3. However, it is relatively low compare d to mean val- 
ues fo und in other cloud core s. For example, Womac k et al.l 
(I1992I) and lBenson et all dl998l) found mean values of X(N 2 H + ) 
of 4 x 10~ 10 and 7 x 10~ 10 , respectively, from various samples. 

From Table [9] we can also determine the average val- 
ues of some relative abundances: X(HCN)/X(HNC) ~ 1.8, 
X(HNC)/X(HCO + ) ~ 0.3, and X(HCN)/X(HCO + ) ~ 0.6. The 
value of Z (HNC)/Z(HCO + ) ag rees reasonable well with that es- 
timated bv lGodard et al.l(l2010l) (0.5+0.3), whereas the other two 
relative abundances, X(HCN)/X(HNC) and X(HCN)/X(HCO + ) 
both appear to be somewhat lower compared to the estimated 
values (4.8 + 1.3 and 1.9 + 0.9, respectively) of Godar d et al.l 
(1201 Ol) . though all estimates have large (<; 50%) uncertainties. 

4.2. Radial profiles of relative molecular abundances 

The computation of the absolute molecular abundances is af- 
fected by our ability to determine the mass distributions of both 
a specific molecular tracer and of the total mass. In particular, 
the spatial distribution (i.e., pixel-by-pixel) of the total mass, as 
determined from the dust continuum, is seriously affected by our 
ability to estimate and subtract the local background emission. 

In order to avoid these difficulties, we have decided to eval- 
uate only the spatial distribution of the relative molecular abun- 
dances, which depend solely on our Mopra observations and are 
not affected by variable levels of background in the bolome- 
ters maps. In Figure [7] we show the source-averaged, radial 
profiles of the normalized relative abundances [HCO + ]/[HCN], 
[HNC]/[HCN] and [N 2 H + ]/[HCO + ] with associated standard- 
deviations. We show separately the resulting averages for proto- 
stellar and starless cores with available map data. For each sepa- 
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Fig. 7. Source-averaged radial profiles of the relative molecular abundances (from left to right) of [HCO + ]/[HCN], [HNC]/[HCN] 
and [N2H + ]/[HCO + ]. The abundances are normalized with respect to the peak value. The top and bottom panels show the resulting 
averages for all starless and proto-stellar cores, respectively, where data were available. Pixel size as in Figure [6j 



Table 10. Abundance estimates. 



Source # 


N 2 H + 


HNC 


HCO + 


HCN 


NH 3 fl 




[xlO- 9 ] 


[xlO- 9 ] 


[xlO" 9 ] 


[xlO- 9 ] 


[xlO- 8 ] 


9 




0.85 


0.50 


0.09 


25.4 


13 




0.57 


1.18 


0.54 




24 


0.06 


0.26 


0.55 


0.18 


0.6 


31 






0.56 


0.16 


96.9 


34 




0.47 


0.41 






40 


0.02 


0.78 


0.22 


0.40 


3.3 


41 






2.63 




2.7 


47 




0.86 


1.03 


2.21 


18.3 


50 






2.59 


4.3 




56 




0.07 


0.84 


0.33 




57 






0.62 






63 


0.11 


0.35 


0.77 


0.62 


6.7 


81 


0.03 


0.19 


0.58 


0.20 


0.08 


82 


0.15 


0.49 


1.8 


0.26 


0.9 


89 






2.80 






101 




0.46 


1.65 


0.44 




109 






0.89 







Notes. The table lists the abundance values [MOL]/[H 2 ] for the main 
observed molecules toward the sources that were mapped. 
(o) Same sources as in Table|7] 



rate source, the relative abundances are normalized with respect 
to the observed maximum value, as a function of radius. We then 
take the median values of these normalized relative abundances 
in each sub-sample of sources. The advantage of using the nor- 
malized values is to be less affected by the variations of the ab- 
solute values of the abundances from source to source, and high- 
lights instead the behaviour of the abundances as a function of 
radius. 

Therefore, in Figure [7] we note that while in the proto- 
stellar cores the [HCO + ]/[HCN] ratio peaks at the center of 
the sources, in the starless cores the peak is slightly off-center. 
This shows that on average the abundance of [HCO + ] relative 
to that of [HCN] is somewhat lower toward the center of the 
starless cores. As far as the [HNC]/[HCN] ratio is concerned, 



our data show that this relative abundance is quite fiat in both 
starless and proto-stellar cores, though a slight decrease as a 
function of radius can be observed, particularly toward proto- 
stellar cores. Finally, we can note the very different behaviour 
of the [N2H + ]/[HCO + ] ratio in starless and proto-stellar cores. 
In fact, while in the former this abundance ratio is quite fiat, in 
proto-stellar cores the [N2H + ]/[HCO + ] ratio can be observed to 
decrease toward larger radii. 

A detailed comparison of these radial profiles with simu- 
lations obtained using re cent chemical models (e.g., iLee et all 
12004 lAikawa et al.l f2008). which follow the chemical evolution 
of a source both as a function of time and of radial position, is 
difficult because they usually model the chemistry of the source 
only out to a radius ~ 1 - 2 x 10 4 AU, a size lower than that 
(26600 AU) covered by the Mopra beam at the distance of Vela- 
D. 

However, we note that some featur es present in Fig ure [7] 
are consistent with, e.g., the model of ILee etal (2004). The 
faster drop-off at large angular radii of the radial profile of the 
[N2H + ]/[HCO + ] ratio in proto-stellar cores is consistent with 
the molecular abundance profiles after collapse discussed by 
ILee et al] d2004l) . though the y are not consistent with similar ra- 
dial profiles discussed by lAikawa et al.l (|2008). Also, the fact 
that in Figure [7] the peak of the [HCO + ]/[HCN] ratio in starless 
cores is slightly off-center could be a consequence of HCO + de- 
pletion, which is more effective in decreasing the column density 
of HCO + toward the c enter of the sourc es and is more effective 
at earlier times. In fact. lLee et al J (12004) have shown that deple- 
tion initiates during the pre-collapse phase and proceeds through 
the collapse of the source. 

4.3. Chemical status of the cores 

In this section we further discuss some chemical implications 
of the relative molecular abundances. It is often found in the 
literature that some molecules are classified either as "early- 
time" (10 4 - 10 6 years after the onset of chemistry) or "late- 
time" species (maximum abundances reached at steady state, af- 
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z 
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Fig. 8. Column density correlations for the same molecules ex- 
amined in Figure |7] The solid line represents the linear fit (from 
Bayesian statistics, see text) to all points. The Spearman rank 
coefficient is > 0.5 in all three cases, with a peak of 0.75 in the 
plot of JV(HNC) vs. JV(HCN). 

ter 10 6 - 10 8 years of chemical evolution), based on the produc- 
tion pathways via ion-molecule or neutral-neutral reactions. Of 
the molecules discussed here HCO + and NH3 are usually consid- 
ered as late-time molecules, but the situation of HCN and HNC 
is less clear in this picture. 

However, we can attempt to use the [HNC]/[HCN] abun- 
dance ratio as a measure of the evolutionary phase and tem- 
perature in our sources. In fact, converting Figure [7] into abso- 
lute values (instead of the normalized values shown in the fig- 
ure) for the relative molecular abundances we find that: ( i) the 
abundances of HCN and HNC are not significantly different be- 
tween the proto-stellar and the starless cores, thus a systematic 
change in the HCN abundance due to chemical evolution is not 
apparent in this study; ( ii) the average [HNC]/[HCN] ratio is in 



agreement with the v alues measured in 19 low-mass starless and 
proto-stellar cores bv lHirota et all(ll998l) . 

According to gas-phase chemical models, HCN and HNC 
are mainly produced by a dissociative r ecombination rea ction 
of the HCNH + ion with an electron (see lHirota etalJI 19981 and 
references therein): 



HCNH + + e 



HCN + H 
HNC + H 



(9) 



The branching production ratio of HCN, a, is defined as 
[HCN]/[HNC]= or/(l - or), if T k <K T c , where T c 24 K is 
the threshold temperature above which neutral -neutral reactions 
dominate the [HCN]/[HNC] ratio. If we take a median value of 
[HCN]/[HNC] - 1 in our cores, the branching production ratio of 
HCN is estima ted to be a ^ 0.5 , somewhat higher than the value 
reported by Hirota et al. ( 199 8J) (a ^ 0.4), but in a greement with 
the range of values found by Nikolic et al.l J2003) toward LI 251 
(a =: 0.2 - 0.8). The relatively low value of a suggests that the 
more energetic paths for destruction of HNC, once the temper- 
ature of the gas exceeds the critical temperature, are less favor- 
able, and thus it constitues a further indirect evidence that our 
cores are cold. 

Another parameter to assess the chemic al status of our cores 
is the [N 2 H + ]/[H 13 CO + ] abundance ratio. iFuente et all d2005h 
have proposed this ratio as a measure of the evolutionary phase 
of a source. According to these authors, since the abundance 
of N 2 H + tends to remain constant in starless clumps, while 
H 13 CO + could suffer from depletion in the densest part, then 
the [N2H + ]/[H 13 CO + ] ratio could be used as an indication of 
the chemical evolutionary phase of the cores. If we use the 
single-point observations of both N 2 H + and H 13 CO + we find 
a median value of their relative abundance, [N 2 H + ]/[H 13 CO + ] = 
5.7 + 1.9, which is inte rmediate in the range of values proposed 
bv IFuente etall (120051) spanning from Class objects (where 
molecular depletion is significant and the [N2H + ]/[H 13 CO + ] ra- 
tio is maximum, ~ 15) to more evolved objects where the ratio 
is minimum (~ 3). Note, however, that the objects studied by 
IFuente et al.l Q005) are on average more evolved than those in 
our sample of cores in Vela-D. 

4.4. Correlations between parameters 
AAA. Column density correlations 

We can use the molecular column densities estimated in the pre- 
vious sections to compare the behavior of the molecules. If we 
analyze molecules with different critical densities, we would ex- 
pect their molecular transitions to be emitted in different vol- 
umes of gas, in the presence of density gradients. Likewise, 
molecules with similar critical densities will be excited in the 
same volume of gas, in the absence of strong abundance gradi- 
ents. 

We performed this test only in those sources where spec- 
tral line maps were available, to ensure that we integrate over 
all molecular emission. Thus, in these correlation plots all in- 
formation about the specific spatial distribution of the different 
molecular tracers is lost. In addition, the inclusion in this analy- 
sis of spectral line maps only, dramatically decreases the size of 
the sample, which also represents a range of masses, chemistries, 
and evolutionary states. 

Our results are presented in Figure[8] where we have used the 
Bayesian IDL routine LINMIX_ERR to perform a linear regres- 
sion, to find the slope of the best fit line in each case. In the top 
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Fig. 9. Median core size (left panels) and median linewidth, AV (right panels) vs. n e s (at 15 K), shown separately for a sub-sample 
of proto-stellar (top row) and starless (bottom row) cores. The plot of AV vs. n e s makes use of the single-point observations, instead 
of spectral line maps as in the left panels, allowing for more molecular lines to be used. The effective densities are listed in Table|2] 
The plot also includes some detections of the HC3N(1 - 0) line, with n e s = 2.5 x 10 4 . The solid line represents the linear fit (from 
Bayesian statistics) to all points. 



panel we plot the column density of HCO + vs. the column den- 
sity of HCN. HCN should be a higher-density gas tracer, com- 
pared to HCO + ; however, despite a considerable scatter, the two 
column densities show a modest correlation, with a Spearman 
rank coefficient > 0.5. The column densities of HCN and its 
isomer HNC also appear to be well correlated (correlation co- 
efficient = 0.75) in the middle panel of Figure [8] Finally, in 
the bottom panel, we also compare the column densities of the 
chemically opposite species N2H + and HCO + . The few points 
in this panel are a consequence of the fewer available maps for 
N2H + compared to the other molecular tracers. However, we still 
see a correlation between the two column densities (correlation 
coefficient > 0.5). 

The variable degree of correlation and the scatter in Figure[8] 
is an indication of source-to-source variations in the physical and 
chemical conditions. However, these qualitative trends are re- 
markable given the diversity of our sources and relatively small 
sample size. 

4.4.2. Correlations with « e ff 

The molecular transitions observed by us are characterized by a 
range of critical densities, thus allowing us to examine how spe- 
cific physical parameters change as a function of the gas traced 
by a specific line. However, instead of using the critical density, 
n cr ; t , in our analysis we prefer to use the effective density, n e s, 



to characterize the volume of gas where a given spectral line is 
generated. As discussed in lEvansI d!999l) . the ability to detect a 
transition is usually associated with a number density larger than 
n cr it, but lines are also easily excited in subthermally populated 
gas with densities more than an order of magnitude lower than 

n crit- 

To account for effects such as optical depth, multile vel ex- 
citatio n effects, and trapping, we use the definition of lEvansI 
(Il999h . where n s s represents the density needed to excite a 
Tr — 1 K line (where Tr is the observed radiation tempera- 
ture) for a given kinetic temperature, as calculated with a non- 
LTE radiative transfer code. We specifically calculated the effec- 
tive densities of our molecular tracers by using the online ver- 
sion of RADEX dvan der Tak et all 120071) . a non-LTE radiative 
transfer cod e, assuming log(AVAv) = 13.5 cm^/lkms -1 ) (see 
lEvansI fl999l) and a kinetic temperature of 15 K for all species 
(see Tabled 

Figure [9] shows the (deconvolved) core diameter, as deter- 
mined from the molecular intensity maps, as a function of n e ff- 
The observed core diameter is estimated using the area enclosed 
within the half -peak intensity contour and then deconvolved us- 
ing the standard Gaussian deconvolution with the Mopra beam. 
We then take the median value of the size (using the selected 
sources listed in Figure [9| and plot it as a function of n e fj in 
Figure [9] where the error bars represent the median absolute de- 
viations. By comparing the values of the measured core diameter 
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Table 11. Median properties of starless and proto-stellar cores. 



Source Type 


Ram 


n 


FWHM" 


T 


*W b 
M vir 


Ptot 
k 


|NH,J c 
[N 2 H+] 


[H 13 CO+] 




[arcsec] 




[kms->] 


[K] 




[xlO 6 cnr 3 K] 




Starless 


27.0 ± 12.0 


0.8 + 0.7 


1.4 ±0.4 


15.1 ± 1.6 


0.5 ± 0.3 


~ 5 


-10* 


5.7 ± 1.7 


BLAST063 


24.0 ± 6.0 rf 


0.6 ± 0.3'' 


1.5 ±0.1 


12.2 


0.3 


~ 10 


-600 


5.7 


Proto-stellar 


21.0 + 5.9 


1.1 ±0.5 


1.6 ±0.6 


16.8 ± 1.8 


0.8 ± 0.4 


~ 19 


~ 60 


6.4 ± 1.9 



Notes. <n> From optically thin tracers only. <6) Af vir estimated assuming a density distribution p oc r 2 . {c} Mean values based on two (starless) 
and three (proto-stellar) sources. (d) Excluding N2H + . 
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Fig. 10. Plot of average linewidth vs. core temperature. The 
BLAST-derived dust temperature has been used, except in 
those sources where 2\ could be determined from NH3 (see 
Section l3~3l l. In estimating the averages, only the linewidths of 
the optically thin molecular tracers (N 2 H + , H 13 CO + and H 13 CN) 
have been considered. The triangles and squares represent star- 
less and proto-stellar cores, respectively. The filling color code 
represents the relative abundance [HNC]/[HCN] (top panel) and 
[N 2 H + ]/[H 13 CO + ] (bottom panel). The solid line represents the 
linear fit (from Bayesian statistics) to all points. 



in starless and proto-stellar cores we note that for the same value 
of « e ff the starless cores have on average a larger size. 

Figure [9] also shows the source linewidth, AV, measured 
from single-point spectra, as a function of n e g. Because in this 
case we are using the single-point observations, we could in- 
clude in the plot almost all molecular transitions listed in Table|2] 
Furthermore, since linewidths may be broadened by optical 
depth effects, for Figure [9] we selected only those source/line 
cases with t < 1. We note that no significant difference in AV 
is observed, for the same value of n e s, between proto-stellar and 



starless cores, but we find that in starless cores AV and n e s are 
correlated. In fact, as for the case of Figure [8] we have used the 
IDL routine LINMIX_ERR to perform a linear regression, and 
we find a Spearman rank coefficient of 0.97 for the starless cores. 
However, this correlation coefficient drops to 0.23 for the proto- 
stellar cores. 



4.4.3. Linewidth-temperature correlation 

Since there is no single measurement that can determine the "ex- 
act" evolutionary stage of each source, as a further test to help 
estimate the status of the cores in our sample we plot in FigurefTUl 
the linewidth vs. the core temperature. In order to avoid possible 
optical thickness effects, the figure only includes the linewidths 
of the (mostly) optically thin molecular tracers, i.e., N2H + and 
the two isotopologues H 13 CO + and H 13 CN. From this plot we 
also had to eliminate those sources with a poor determination of 
the linewidths because of weak lines and/or noisy spectra. 

Figure[10]thus shows that there is a clear segregation in tem- 
perature between starless and proto-stellar cores. The figure also 
shows that there is a modest correlation between linewidth and 
temperature. As for the case of Figure [8] we have used the IDL 
routine LINMIX_ERR to perform a linear regression, and we 
find a Spearman rank coefficient of 0.47. This modest correlation 
is confirmed by the fact that the median values of the linewidth 
determined for the starless and proto-stellar cores are 1.4 + 0.4 
and 1.6 + 0.6 km s" 1 , respectively (see Table [TTli. which are con- 
sistent within the errors. 

Figure [10] also encodes the information about the relative 
molecular abundances [HNC]/[HCN] and [N 2 H + ]/[H 13 CO + ], 
discussed in Section |3] No apparent trend of [HNC]/[HCN] 
or [N2H + ]/[H 13 CO + ] with either linewidth or temperature can 
be seen. 



4.5. Comments to specific sources mapped with the Mopra 
telescope 

In this section we further discuss some starless and proto-stellar 
cores that have been mapped at Mopra, including the "transition" 
source BLAST063. 

BLAST024. This source is a compact proto-stellar core, 
where the spatial distribution of the emission of all main trac- 
ers observed by us is centered on the BLAST dust core. Wing 
emission is detected in the single-point HCO + (l - 0) spectrum, 
suggesting the presence of a molecular outflow. However, due 
to the low SNR in our maps, we could not produce an image 
of the outflow, as it is the case also for sources BLAST08 1 and 
BLAST082. 

BLAST031. This is a moderately compact starless core, 
where the HCN(1 -0) and HCO + (l -0) emission trace quite well 
the BLAST continuum emission, but the HNC( 1-0) emission is 
seen to peak at an offset position compared to the BLAST emis- 
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sion and the other two molecular probes. Such offsets are also 
found toward other sources (see below). 

BLAST034. This is an irregularly shaped, starless core 
which shows two well separated velocity components in the 
single-point HCO + (l - 0) spectrum, with the main (mapped) 
component being the low-velocity one (^ 2kms~ 1 ). The two 
velocity components are also detected in the HNC(1 - 0) spec- 
trum. 

BLAST047. This source also shows a double-peaked 
HCO + (l - 0) spectrum. However, both the channel map of the 
HCO + (l - 0) emission and the spatial distribution of the HNC 
and HCN(1 - 0) emission (see Figure lA~2l suggest that HCO + is 
self-absorbed at the position of the BLAST peak. One can also 
note that HCN, and especially HNC, trace very well the BLAST 
contour plots. 

BLAST063. This source is the most intense starless core 
mapped by us at Mopra. This core has been classified throughout 
this pa per and was originally classified as starless bv lOlmi et al.l 
(2009) because no MIPS compact source at 24 \mi could be 
found at its location. However, we have found that the single- 
point HCO + (l - 0) spectrum shows line-wing emission, indi- 
cating the presence of a molecular outflow, which should actu- 
ally change the source type from starless to proto-stellar. The 
fact that no 24 ymi emission is detected toward this core, though 
a molecular outflow is already active, suggests that the proto- 
stellar core has just formed and is thus in an earlier stage com- 
pared to BLAST024, BLAST081 and BLAST082. Therefore, 
this core could tentatively be classified as being in a "transition 
phase" from starless to proto-stellar. We also note that this source 
has an elongated shape, and we can see that the HCN(1 - 0) 
and HNC(1 -0) emission trace very well the BLAST continuum 
emission, whereas there is an offset between the HCO + (l - 0) 
emission and the dust emission. 

BLAST081. This source is another compact proto-stellar 
core mapped with the Mopra telescope. Despite some artefacts 
visible in the HCO + (l - 0) map (see Figure lA~4i i we note that 
while the N 2 H + (1 - 0) and (possibly) HCO + (l - 0) emission 
follow closely the BLAST continuum emission, the emission of 
the two isomers HCN and HNC peak at an (identical) offset po- 
sition compared to the BLAST maximum. This is the only ex- 
ample of such an offset observed among the three proto-stellar 
cores mapped with the Mopra telescope, though it is also partly 
observed toward the early proto-stellar core BLAST063. Since 
our sensitivity was not high enough to map the optically thin 
isotopologues, such as H 13 CO + and H 13 CN, we cannot deter- 
mine at present whether this offset is due to optical depth ef- 
fects or chemical variations. As in BLAST024, the HCO + (l - 0) 
spectrum shows line-wing emission, indicating the presence of a 
molecular outflow. 

BLAST082. This source is the third compact proto-stellar 
core mapped at Mopra. Contrary to BLAST081, all mapped 
molecular tracers are well centered on the BLAST continuum 
emission and also follows quite well the less dense material. 
However, the HCN(1 - 0) and HNC(1 - 0) maps show an ad- 
ditional smaller core, located at (Al, Ab) (-60", -60") (see 
Figure IA.4t . which is not visible in the HCO + (l - 0) map. 
This sub-structure may be a consequence of core fragmenta- 
tion and the difference observed in the various molecular trac- 
ers may indicate a different chemical status compared to the 
main core. Then, like BLAST034 and BLAST047, the single- 
point HCO + (l -0) spectrum is also double-peaked, but the chan- 
nel map suggests this is the consequence of two velocity com- 
ponents aligned along the line of sight. Like BLAST024 and 
BLAST08 1 the HCO + ( 1 -0) spectrum also suggests the presence 



of a molecular outflow. Line wings in the HCO + (l - 0) spectrum 
are also observed toward other sources, as listed in Table IC. 41 

4.6. Median properties of starless and proto-stellar sources 
and comparison with other surveys 

In this section we summarize the main differences between star- 
less and proto-stellar cores, by presenting a list of their me- 
dian properties in Table QT| For comparison, we have also in- 
cluded source BLAST063 which we have defined as a transition 
source between the two classes of objects. One can note that 
the starless cores appear to be moderately colder and less turbu- 
lent compared to the proto-stellar sources, though the differences 
are within the uncertainties. The linewidths have been averaged 
among all optically thin tracers (N 2 H + (1 - 0), H 13 CO + (l - 0) 
and H 13 CN + (1 - 0)) and all sources. The lower temperature of 
the starless cores is simply a consequence of the analysis car- 
ried out by lOlmi et aTl (2009), but we now find that also their 
linewidths do not strongly differ from those of the proto-stellar 
cores. 

We do find a more significant difference between starless and 
proto-stellar cores when we consider the radial profiles of the 
column densities. By using the results of Section [3~6l the median 
values of the parameters /?a at confirm that the proto-stellar cores 
are more compact on average. The parameter 77, on the other 
hand, is subject to larger uncertainties, also due to our relatively 
low spatial resolution at the distance of Vela-D, and no definitive 
statement can thus be made. We note that the larger Rf\. dt mea- 
su red toward th e starle ss cores agrees well with the conclusion 
of lCaselli et all (120021) . In fact, these authors found that while 
proto-stellar cores in their sample (which includes sources from 
various star forming regions) were better modeled with single 
power-law density profiles, the starless cores presented a central 
flattening in the integrated intensity profile. 

In terms of the kinem atics of the cores, ou r findings are also 
very similar to those of ICaselli et al.1 (120021) . In fact, we find 
an average value for the velocity gradient in our sample (from 
the HCO + results listed in Table |6} of 1.4 + 0.3 kms -1 pc _1 , 
and we do not find significant differences between starless and 
proto-stellar cores. This is consistent with t he typical value of 
2+1 kms 1 pc~' found bv lCaselli et al.l(l2002|) in their sample. In 
terms of the equilibrium status of the cores, ICaselli et al ] (l2002l) 
found that their "excitation mass" to virial mass ratio is typically 
M ex /M w \ r si 1.3, for both classes of cores. From Table [TT1 on the 
other hand, we can see that both starless and proto-stellar cores 
have 0.5 ^ M core /M vil ^ 1. Therefore, while the M mre /M v i r ra- 
tio is somewhat lower compared to the results of lCaselli et al.l 
(2002) (but this may also be a consequence of the different meth- 
ods to calculate the core masses) it is still very near or beyond 
the self-gravitating threshold of 0.5. We also find that in proto- 
stellar cores the M core /M v ; r ratio is slightly larger than in starless 
cores, though they are still consistent within the errors. 

In Table QT] we have also included the median values of the 
average total internal pressure (not the edge pressure), P tot , of 
each core in Vela-D for which we have measured (optically thin) 
linewidths. The values of P tr ,t have been ca lculated following 
the method described by lOlmi et aTl (2010) and are given in 
the usual cm 4 K units for P to i/k. We note that in proto-stellar 
cores the median internal pressure is much higher than in star- 
less cores, as expected. However, the confinement of the cores 
is of concern only for those sources (mostly starless) that are 
gravitationally unbound. Therefore, the starless cores which are 
gravitationally unbound according to the M core /M v ; r ratio can be 
stable (against expansion) only if their edge pressure is less than 
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the local a mbient pressure (o f order P ext /k ~ 5 x 10 5 cm 3 K in 
Vela-D, see lOlmi et alJ20loh . Pressure confined cores have been 
previously discus sed, for example, by lLada et all d2008l) and 
ISaito et alJ d2008l) . Alternativ e sources of supp ort could come 
from the magnetic field (see lOlmi et al] 1201 Ol and references 
therein), or the gravitationally unbound cores must otherwise be 
transient structures. 

As far as chemical properties are concerned, TablefTTIshows 
that the values of the relative abundance [N 2 H + ]/[H 13 CO + ] for 
the starless and proto-stellar cores are quite similar. Thus, as 
noted in Section 1431 we do no t find any trend bas ed on this ra- 
tio, as it has been suggested bv lFuente et al.l d2005l) . However, in 
Table [TT1 we note that the [NH3]/[N2H + ] abundance ratio is sig- 
nificantly larger toward starless cores, though we had only two 
starless cores where we could effectively measure this ratio. This 
trend, if not the quantitative difference between the two classes 
of objects ( because of th e low s tatistics), is consistent with the 
findings of iFriesen et al.l d2010l) . who concluded that the rela- 
tive fractional abundance of NH3 to N2H + remains larger toward 
starless cores than toward protostellar cores by a factor of ~ 2-6. 

These properties, and in particular the fact that many of the 
starless cores are gravitationally bound (hence, pre-stellar), sug- 
gest that at least most of the pre-stellar cores observed by us 
are indeed the precursors of the proto-stellar cores and not just 
transient objects following a different evolutionary path. Their 
similar median temperature and linewidths indicate that the en- 
velopes of the sources have not had enough time to "register" the 
emergence of a protostar inside the core. In turn, this also means 
that most of the turbulence in the proto-stellar cores observed by 
us is generated prior to the appearance of a warm object inside 
the core. 

Therefore, it appears that the transition from the pre- to the 
proto-stellar phase is relatively fast, leaving the core envelopes 
with almost unchanged physical parameters. Alternatively, if 
most proto-stellar cores observed had not had much time to 
evolve after the appearance of a protostar, this might also explain 
the similar median temperatures and linewidths of the two source 
types. However, the latter hypothesis seems less likely given that 
among our s elected proto-stel lar cores are also sources relatively 
evolved (see Olmi et al. 2009 and references therein). 

In this respect, despite the uncertainties it is of interest to 
note that the values of the parameters /?g at , 77 and FWHM of 
source BLAST063 are somewhat intermediate between those of 
starless and proto-stellar cores. But, if the transition phase from 
pre- to proto-stellar is indeed fast, then one should not expect to 
detect many such transition objects, and their properties may be 
indistinguishable from those of either class of sources. More sys- 
tematic observations of larger samples of sources are required to 
further study this issue. 

5. Summary and conclusions 

We presented 3-mm molecular line observations of the Vela-D 
region obtained with the 22-m Mopra radio telescope and 1.3 cm 
NH3 observations carried out with the Parkes antenna. In total 8 
molecular lines were analyzed in 40 sources, both starless and 
proto-stellar cores, previously detected by BLAST. A total of 20 
spectral line maps were obtained. The results of our study can be 
summarized as follows: 

1 . Our spectral line maps show a wide variety of morphologi- 
cal types: very early and cold starless cores appear to have 
an irregular shape (e.g., BLAST009) in most or all molec- 
ular tracers mapped at Mopra. Some warmer cores (e.g., 



BLAST031) and cores at the transition phase from starless 
to proto-stellar (BLAST063) appear to be more regularly 
shaped and more compact. Finally, proto-stellar cores all 
show a more rounded shape and narrow radial intensity pro- 
files. 

2. We compared the kinetic temperature, 7\, derived from the 
NH 3 (1,1) and (2,2) observations, with the BLAST-derived 
dust temperature, Tj, and found that the median value of 
their ratio is 7\/7d =1.2 + 0.2. 

3. We found that all proto-stellar sources, and at least 7 starless 
cores mapped at Mopra, show velocity gradients. In two of 
the proto-stellar cores the direction of the velocity gradients 
measured in HNC and HCO + are significantly different, sug- 
gesting the possibility that different systematic motions are 
simultaneously present (one of which likely is a molecular 
outflow). 

4. The analysis of the virial masses showed that nearly all of the 
starless cores have masses below the self-gravitating thresh- 
old, indicating that they are unlikely to be gravitationally 
bound, whereas more than half of the proto-stellar cores have 
masses which are near or above the self-gravitating critical 
value. If a density profile closer to the one actually observed 
is considered in the calculation of the virial masses, then also 
most starless cores turn out to have masses near or above the 
self-gravitating threshold. 

5. The average internal pressure of proto-stellar cores is higher 
than in starless cores, and both are higher than the typical 
ambient pressure in Vela-D. Cores (mostly starless) not grav- 
itationally bound cannot thus be pressure-confined like, for 
example, the cores in the Pipe nebula. 

6. The radial profile of the ring-averaged integrated spectral 
line intensity confirmed that proto-stellar cores have on av- 
erage a more compact structure. We also found that the 
radial profile of the N2H + (1 - 0) emission falls-off more 
quickly, on average, than that of the C-bearing molecular 
lines HNC(1 - 0), HCO + (l - 0) and HCN(1 - 0). 

7. We find a variable degree of correlation, with a significant 
scatter, between the column densities of chemically differ- 
ent molecular species, and also between molecular tracers 
with different effective densities. We also find that in star- 
less cores the linewidth and effective density are well cor- 
related. Instead, linewidth and temperature are not strongly 
correlated. 

8. The branching production ratio of HCN and also the rela- 
tively low [N2H + ]/[H 13 CO + ] abundance ratio confirm that 
all cores in our sample are cold and are on average in early 
evolutionary phases. 

9. An analysis of the median properties of the starless and 
proto-stellar cores suggests that the transition from the pre- 
to the proto-stellar phase is relatively fast, leaving the core 
envelopes with almost unchanged physical parameters. 
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Fig. A.l. Maps of the line integrated intensities (from left to right, N 2 H + (1 - 0), HCN(1 - 0), HCO + (l - 0) and HNC(1 - 0), in 
units of J T* dv [Kkms 1 ]) of selected sources toward Vela-D. In this figure we show, from top to bottom, maps of BLAST009, 
BLAST013, BLAST024 and BLAST031. The white dashed contours represent the BLAST flux density at 250yum. 
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Fig. A.2. Same as Figure|AT]for sources BLAST034, BLAST040, BLAST041 and BLAST047. 
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Fig. B.l. Spectra of N 2 H + (1 - 0) towards the Vela-D sources. 
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B.2. Same as Figure|Ofor the HCN(1 - 0) line 



Please give a shorter version with: \authorrunning and/or \titilerunning prior to \maketitle 



BLAST054: HCN(1 -0) 



BLAST056: HCN(1 - 




BLAST057: HCN(1-0) 



BLAST093: HCN(1 -0) 



BLAST097: HCN(1-0) 



2.0 










1.5 






1.0 






1.0 






t 0.5 




\ ■ 


0.5 
0.0 






0.0 







-20-10 
Velocity 



10 20 30 
[km/s] 



-30 -20 -10 10 20 30 40 
Velocity [km/s] 




-30 -20 -10 
Velocity 



10 20 30 40 
[km/s] 



BLAST065: HCN(1-0) 




-30 -20 -10 
Velocity 



10 20 30 40 
[km/s] 



BLAST077: HCN(1-0) 





-30 -20 -10 

Velocity 



10 20 30 40 
[km/s] 



BLAST082: HCN(1-0) 



2.0 - 
, 1.5 - 
' 1.0 r 

0.5 : 

0.0 ^Aw^wyA/V 

-30 -20 -10 
Velocity 



10 20 30 40 
[km/s] 



BLAST090: HCN(1-0) 



-30 -20 -10 10 20 30 40 
Velocity [km/s] 



B.3. Ctnd. from Figure EHl 



Please give a shorter version with: \authorrunning and/or \titilerunning prior to \maketitle 



BLAST003: HNC(1 -0) 



BLAST009: HNC(1 -0) 



BLAST023: HNC(1-0) 





0.4 






0.6 








S 0.2 






„ 0.4 














{ 0.2 








0.0 






0.0 







-30 -20 -10 10 20 30 40 
Velocity [km/s] 



BLAST024: HNC(1 -0) 



— 3 

^ f 2 



o A^yyA**^ 




-30 -20 -10 10 20 30 40 
Velocity [km/s] 

BLAST034: HNC(1 -0) 




-30 -20 -10 10 20 30 40 
Velocity [km/s] 



BLAST041 : HNC(1 -0) 




-30 -20 -10 10 20 30 40 
Velocity [km/s] 



BLAST045: HNC(1 -0) 




-30 -20 -10 10 20 30 40 
Velocity [km/s] 



BLAST052: HNC(1 -0) 




-30 -20 -10 10 20 30 40 
Velocity [km/s] 



-30 -20 -10 10 20 30 40 
Velocity [km/s] 

BLAST026: HNC(1-0) 



— o.i 




-30 -20 -10 10 20 30 40 
Velocity [km/s] 

BLAST038: HNC(1 -0) 



0.3 
0.2 
0.1 
0.0 
-0.1 



-30 -20 -10 10 20 30 40 
Velocity [km/s] 

BLAST043: HNC(1-0) 



0.6 

0.4 

S 0.2 
I 

0.0 

-0.2 




-30 -20 -10 10 20 30 40 
Velocity [km/s] 

BLAST047: HNC(1 -0) 



0.3 
0.2 
0.1 
0.0 
-0.1 




-30 -20 -10 10 20 30 40 
Velocity [km/s] 

BLAST053: HNC(1 -0) 



0.6 
' 0.4 
0.2 
0.0 



-30 -20 -10 10 20 30 40 
Velocity [km/s] 



-30 -20 -10 10 20 30 40 
Velocity [km/s] 

BLAST031: HNC(1-0) 





„ 
1 



-0 



-30 -20 -10 10 20 30 40 
Velocity [km/s] 

BLAST040: HNC(1-0) 



Mi 



-30 -20 -10 10 20 30 40 
Velocity [km/s] 

BLAST044: HNC(1-0) 



0.0 'kuWfSjfil 



-30 -20 -10 10 20 
Velocity [km/s] 



BLAST050: HNC(1-C 




-30 -20 -10 10 20 
Velocity [km/s] 



BLAST054: HNC(1~C 




-30 -20 -10 10 20 
Velocity [km/s] 
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Fig. B.6. Same as Figure |BTT] for the HCO + (l - 0) (lower spectra) and H 13 CO + (l - 0) (upper spectra) lines. Only sources where 
H 13 CO + (l - 0) was detected are shown. 
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Fig. B.7. Same as Figure iBTl for the H 13 CN(1 - 0) (upper spectra) lines, shown together with the corresponding HCN(1 - 0) (lower 
spectra) lines. Only sources where H 13 CN(1 - 0) was detected are shown. 
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Fig. B.8. Spectra of the NH 3 (1,1) line obtained with the Parkes telescope. 
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Appendix C: Molecular lines parameters 
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Table C.l. N2H + line parameters. 



Source # 






N 2 H + 










Vlsr 


FWHM 


T 


T„ 




[K] 


[km s" 1 ] 


[kins" 1 ] 




[K] 


3 


0.57 


4.3 


1.0 


0.72 


4.5 


9 


0.35 


9.8 


0.7 


0.10 


- 


24 


2.84 


3.2 


2.3 


0.10 


- 


40 


0.79 


0.3 


0.7 


0.10 


- 


44 


0.79 


0.2 


0.9 


0.25 




45 


0.28 


10.4 


1.1 


0.10 




63 


0.34 


11.2 


1.6 


0.78 


3.7 


81 


3.05 


11.6 


2.7 


0.10 




82 


1.38 


10.8 


1.6 


0.65 


7.2 


90 


1.51 


10.4 


1.4 


0.90 


6.3 


93 


1.40 


1.5 


1.5 


0.72 


6.9 


97 


0.55 


7.7 


1.6 


1.09 


3.8 



Notes. Only sources with a single-point observation are listed (see TablefT}. Undetected sources and unreliable hfs fits are not shown. 
Table C.2. HCN and H 13 CN line parameters. 
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[Kill S J 




3 


0.83 


7.1 


2.8 


0.12 












24 


3.12 


5.1 


3.6 


0.16 




0.25 


2.8 


2.9 


0.18 


26 


0.06 


3.1 


1.7 


0.10 












26 


0.06 


3.1 


1.7 


0.10 
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0.22 


11.4 


2.0 


2.06 


3.0 










38 


1.17 


3.6 


1.2 


17.03 


2.9 










40 


0.43 


2.2 


1.3 


0.83 


3.8 










41 


0.29 


12.1 


1.4 


0.10 












43 


0.21 


9.7 


5.2 


1.77 


3.0 










44 


0.49 


2.1 


1.6 


0.84 


4.0 










45 


0.24 


12.1 


1.9 


0.10 












47 


0.50 


6.4 


1.8 


3.08 


3.1 










50 


0.93 


7.7 


2.9 


13.81 


2.9 










52 


0.36 


3.3 


1.1 


4.13 


2.9 










53 


0.22 


10.7 


2.6 


0.10 












54 


0.42 


5.4 


0.9 


3.21 


3.0 










56 


0.16 


10.4 


1.9 


0.10 












57 


0.38 


14.2 


1.5 


3.06 


3.0 










59 


0.23 


8.4 


1.9 


0.53 


3.7 










63 


0.78 


12.7 


2.4 


1.11 


4.2 










77 


0.28 


2.9 


1.8 


0.10 












79 


0.36 


12.1 


2.1 


2.02 


3.1 










81 


4.05 


13.2 


3.1 


0.10 




0.38 


11.3 


3.1 


0.10 


82 


0.67 


12.8 


3.0 


0.11 












88 


0.27 


23.7 


1.8 


0.10 












89 


0.35 


3.1 


1.5 


2.08 


3.1 










90 


1.27 


12.1 


2.1 


0.10 




0.12 


10.2 


1.6 


0.18 


93 


0.80 


3.1 


2.2 


0.24 




0.11 


1.1 


1.5 


0.10 


97 


0.53 


9.1 


2.7 


0.10 













Notes. Only sources with a single-point observation are listed (see Table fTJ. Undetected sources and unreliable hfs fits are not shown. If two 
velocity components were detected in a given source, only the more intense component is listed here. 
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Table C.3. HNC line parameters. 



Source # 




HNC 








v hr 


FWHM 




[Kkms- 1 ] 


[kms 4 ] 


[kms- 1 ] 


3 


1.79 


6.9 


2.9 


9 


0.31 


12.2 


1.5 


23 


0.42 


2.7 


1.6 


24 


7.68 


5.5 


3.2 


26 


0.12 


3.1 


1.4 


26 


0.08 


14.7 


0.8 


31 


0.25 


11.7 


1.6 


34 


0.21 


2.7 


2.9 


38 


0.21 


3.4 


1.3 


40 


0.66 


2.4 


1.4 


41 


0.23 


12.1 


1.1 


43 


0.42 


11.1 


1.5 


44 


0.90 


2.4 


1.5 


45 


0.44 


12.6 


1.7 


52 


0.20 


3.7 


1.9 


53 


0.53 


11.3 


1.8 


54 


0.22 


5.7 


1.7 


56 


0.14 


10.6 


1.6 


57 


0.61 


14.7 


2.8 


59 


0.42 


8.6 


1.8 


63 


1.38 


13.2 


2.4 


72 


0.36 


12.8 


2.4 


77 


0.46 


3.3 


1.8 


79 


0.55 


12.7 


2.3 


81 


8.07 


13.6 


3.1 


82 


1.97 


13.1 


2.5 


88 


0.30 


24.0 


1.5 


89 


0.13 


3.6 


1.4 


90 


2.71 


12.6 


2.0 


93 


2.01 


3.7 


1.8 


97 


1.64 


9.7 


2.5 



Notes. Only sources with a single-point observation are listed (see TablefTJ. Undetected sources and unreliable fits are not shown. If two velocity 
components were detected in a given source, only the more intense component is listed here. 
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Table C.4. HCO + and H 13 CO + Line Parameters. 



Source # HCCF H 1J CO + 





nr 


v br 


FWHM 


T 


T„ 


Line Wings? 






FWHM 


T 




[K] 


[kms 1 ] 


[kms- 1 ] 




[K] 




[K] 


[kms" 1 ] 


[kms" 1 ] 




3 


2.66 


6.3 


3.6 


7.97 


4.4 


NG 


0.33 


5.6 


2.4 


0.20 


9 


0.45 


11.1 


1.9 


12.07 


3.3 


NG 


0.17 


11.3 


2.8 


0.30 


23 


0.37 


1.2 


0.9 


- 


3.7 


N 


— 


- 


- 


— 


24 


13.00 


4.4 


3.8 


4.13 


10.1 


Y 


0.99 


4.7 


2.9 


0.10 


26 


0.31 


1.9 


2.0 


- 


3.1 


N 


- 


- 


- 


- 


26 


0.20 


13.5 


1.5 


- 


3.1 


N 


- 


- 


- 


- 


31 


0.60 


10.4 


1.7 


- 


3.6 


N 


- 


- 


- 


- 


34 


0.50 


2.0 


2.0 


- 


3.3 


N 


- 


- 


- 


- 


34 


0.24 


13.1 


1.6 


- 


3.1 


N 


- 


- 


- 


- 


38 


0.36 


1.9 


1.8 


- 


3.2 


N 


- 


- 


- 


- 


40 


1.30 


1.3 


1.7 


- 


4.5 


N 


0.23 


1.3 


0.8 


- 


41 


0.62 


11.0 


1.4 


- 


3.7 


N 


- 


- 


- 


- 


43 


0.63 


9.9 


1.4 


- 


3.8 


N 


- 


- 


- 


- 


44 


1.85 


1.2 


2.0 


6.65 


4.8 


N 


0.16 


1.5 


1.1 


0.17 


45 


1.14 


11.1 


2.1 


- 


4.0 


NG 


- 


- 


- 


- 


47 


0.37 


5.4 


2.5 


- 


3.1 


N 


- 


- 


- 


- 


47 


0.18 


9.2 


1.3 


- 


3.1 


N 


- 


- 


- 


- 


52 


0.44 


2.3 


2.3 


- 


3.2 


N 


- 


- 


- 


- 


54 


0.24 


4.7 


1.9 


- 


3.1 


N 


— 


- 


- 


- 


56 


0.35 


9.4 


2.0 


- 


3.2 


N 


- 


- 


- 


- 


57 


1.04 


13.2 


2.9 


- 


3.6 


N 


0.21 


14.1 


1.5 


- 


59 


0.51 


7.3 


2.1 


- 


3.3 


N 


0.08 


7.8 


0.9 


- 


63 


2.04 


11.7 


2.5 


4.99 


4.6 


Y 


0.13 


11.8 


1.4 


0.12 


71 


0.17 


7.5 


2.1 


- 


3.0 


N 


- 


— 


- 


- 


72 


0.86 


1.9 


2.1 


- 


3.7 


N 


- 


- 


- 


- 


72 


0.58 


11.6 


2.9 


- 


3.2 


N 


- 


- 


- 


- 


77 


1.00 


1.8 


2.1 




3.9 


N 










79 


0.87 


11.0 


3.8 




3.3 


N 










81 


16.33 


12.4 


3.6 


2.07 


13.7 


Y 


0.95 


12.5 


3.6 


0.05 


82 


0.66 


10.0 


1.1 




4.1 


Y 


0.44 


11.9 


1.8 




82 


2.52 


12.6 


2.5 


11.49 


5.0 


Y 


0.44 


11.9 


1.8 


0.29 


89 


0.57 


2.3 


2.4 




3.3 


N 










90 


5.37 


11.4 


2.3 


7.57 


7.8 


Y 


0.65 


11.5 


1.6 


0.19 


93 


3.52 


2.1 


2.3 


14.82 


6.1 


Y 


0.73 


2.5 


1.5 


0.37 


97 


1.81 


8.1 


2.2 


6.64 


4.6 


Y 


0.28 


8.7 


2.2 


0.17 



Notes. Only sources with a single-point observation are listed (see TablefT). Undetected sources and unreliable hfs fits are not shown. In the 7' 
column "NG" indicates a non-Gaussian line profile. 
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